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Abstract 

We present model equations for the trapping and accumulation of particles in a cylindrical channel with 
nanostructured inner walls when a fluid passes through, carrying a moderate load of impurities. The basic ingredient 
of the model is the introduction of a phenomenological 'effective-charge density' of the walls, related to the electrical 
charges exposed in the nanotexture. The effective charge is gradually reduced as the flow runs through the channel 
and the trapped impurities cover the internal walls. Based on the proposed equations, the position and time 
dependence of the areal density of trapped impurities, and the filtration performance, may be calculated. It is 
proposed that experimentally testing these results may help to understand the enhanced trapping capability 
observed in many diverse nanotextured channel structures. 

Keywords: Nanostructures embedded in larger systems, Fluid impurity trapping in nanostructures, Effective-charge 
density 



Background 

Recently, [1-8] researchers from both academia and indus- 
try have experimentally demonstrated that a variety of 
nanostructures and nanotextured media (see some exam- 
ples later in this introduction) can efficaciously trap 
nanoimpurities carried by fluids when a flow is induced 
by external hydrostatic pressure. These findings are not 
only scientifically interesting, but also promising for the 
socially and economically important application of purifi- 
cation of drinking water and other liquids [4,7-9]. When 
compared to conventional porous niters, the new media 
have the important advantages of retaining impurities 
of sizes typically in the tens of nanometers and, at the 
same time, presenting a resistance to hydrodynamic flow 
orders of magnitude smaller than what conventional mod- 
els would predict for channels of diameters as small as the 
particles being trapped. 

Roughly, we can divide the structures presenting such 
enhanced impurity trapping capability into two groups: 
(a) The first group corresponds to those formed by 
nanometric-diameter channels through which the fluid 
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flows [1-4]. A well-known example is the nanotube 
arrays grown and experimentally tested by Srivastava 
and coworkers [1]. Other specially interesting examples 
are graphene membranes although, by now, they have 
been probed only through molecular dynamics simula- 
tions [2]. In any nanometric-diameter channel, simple 
size exclusion will play a major role in the retention of 
nanoimpurities. However, in addition, these structures 
also exhibit remarkable capability to trap some ions sig- 
nificantly smaller than the channels' diameter [1,2]. The 
resistance to flow is observed to be well lower than what 
conventional models predict for these diameters, a phe- 
nomena often attributed to water-nanostructure interac- 
tions (see, e.g., [1]) though not yet fully understood at the 
quantitative calculation level, (b) The second group corre- 
sponds to nanostructures embedded in larger structures, 
resulting in filters composed by channels with micromet- 
ric diameters and inner walls coated with nanoparticles. 
Examples are conventional microfilters coated with Y2O3 
[5], Zr02 [6], or AI2O3 [7,8] nanopowders (further exam- 
ples can be found in the reviews [3,4,9]). These structures 
have been observed by their growers to have a surprisingly 
good nitration performance for nanometric impurities, as 
small as approximately 10 nm, in spite of the relatively 
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large diameter of the channels (note that in a channel 
with a diameter of 1 /xm only about 0.04% of the fluid 
will transit closer than 10 nm from the walls) [3-9]. Their 
hydrodynamic resistance is quite low, similar to the one 
of conventional micrometric filters. Their trapping capa- 
bility is observed to depend on pH and zeta potential 
[5-8] and, thus, electrostatic and polar attraction may 
be suspected to play a significant role in the filtration 
mechanism and dynamics. However, attempts to modelize 
them have been scarce. The authors of [7,8] empirically 
characterized their niters using general-purpose plug-flow 
adsorption models, like those used for column chro- 
matography, and fitting the Langmuir and BET isotherms. 
While these models may be very useful for the design of 
improved filters, further theoretical work modelling the 
trapping capability and time evolution explicitly in terms 
of the charge of the walls, and also more detailed mea- 
surements, seems to be needed to further understand the 
impurity trapping mechanism and dynamics of this type 
of channels with nanostructured walls. 

The purpose of this 'Nano Idea Letter' is to propose a 
specific model for the nanoimpurity trapping capability of 
cylindrical-like channels with nanostructured inner walls 
of the type composing filters of category 'b' in the previ- 
ous paragraph. We explore theoretically a simplified but 
realistic view in which the improved filtration capability is 
primarily due to the fact that the nanotexturing exposes 
electrical charges in the walls which induce both electro- 
static and van der Waals attractions over the impurities 
in the fluid. This nanostructuring also provides chemical 
anchors for the binding of those impurities once they col- 
lide with the channel walls. Correspondingly, our basic 
ingredients will be the introduction of an effective -charge 
density, z e , of the inner walls of the channels and writing 
down as a function of z e the impurity trapping probability. 
As it could be expected, z e will depend on the areal den- 
sity n of impurities already trapped in the inner walls of 
the channel. We obtain within the model the evolution of 
n and z e with position x and with time t when the liquid is 
flowing through the channel. 

The model produces agreement with the initial trapping 
performances quantitatively reported by experimentalists 
in various systems. Also, we propose that further detailed 
measurements as a function of time may be crucial to 
test these ideas more thoroughly. We believe that some 
aspects of the model could also be useful to partly explain 
the trapping of the smaller ions in the nanodiameter chan- 
nels of category a\ However, its full applicability to that 
case is limited by our use of classical dynamics for the 
carrying fluid. Hence, we do not focus here on that cat- 
egory (also, for these nanodiameter channels, in which 
the number of fluid atoms is manageably small, molecu- 
lar dynamics simulations as those in [2] could be a more 
reliable, albeit not general, approach). 



Obtainment of an equation for the areal density of 
trapped impurities in a channel with 
nanostructured walls 

Initial modelling and notations 

Our starting point, and most of our basic notations, is 
illustrated in Figure 1. We consider a channel with nanos- 
tructured inner walls, its nominal shape being cylindrical- 
like with average radius rn and length L. We divide it into 
slices along the axial coordinate x, each with differential 
thickness dx. A fluid flows through the channel due to 
externally applied hydrostatic pressure, carrying a load of 
impurities. Some of those impurities will become trapped 
by the inner wall of the channel, then reducing its effec- 
tive radius to a value r e (x,t) = rn — r\n{x, t), where 
t is the time, n is the number areal density of trapped 
impurities in the inner wall, and r\ is a constant propor- 
tional to the average volume of impurities. Throughout 
this letter, by areal density' we refer to quantities nor- 
malized using the nominal area of the inner wall ^lizr^olx 
for a differential slice) and not the cross section of the 
channel. Also, for simplicity, we consider all impurities 
equal among them (subsequent generalization to multiple 
chemical species should be easy). The average radius of 
the impurities is noted po. The impurity concentration in 
the fluid is considered to be moderate enough as to not 
significantly affect its viscosity and as for the impurities in 
the fluid to be noninteracting with each other (specially 
when colliding with the channel wall). 




channel nontrapped impurity 



Figure 1 Representation of a nanostructured channel filter as 
modelled in the present letter. The nominal shape of the channel is 
supposed to be cylindrical with length L, and the figure shows only 
the differential slice with axial coordinate from x to x+dx. The radiuses 
r 0 and po correspond to the average dimensions of the bare channel 
and impurities. The effective radiuses r e and p e vary as trapped 
impurities cover the inner wall, via their dependences on, respectively, 
the areal density n of trapped impurities and on the areal density z e of 
effective charge of the inner wall. Thisz e reflects that exposed charges 
in a nanostructured surface attract the impurities in the fluid and also 
constitute binding anchors for those impurities. It is expected to 
diminish as impurities cover the surface, for which we assume the 
simple z e (n) dependence given by Equation 1 of the main text. 
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Effective-charge density of the inner wall, z e 

We now introduce the important concept of a phe- 
nomenological effective charge' of the inner wall of the 
channel We quantify this effective charge via its areal den- 
sity z e , and as already commented on in the introduction, 
it reflects the fact that the nanostructured walls expose 
charges that induce both electrostatic and van der Waals 
attractions over the components of the impurities in the 
fluid. Indeed, z e will depend on the areal density of already 
trapped impurities n (which will screen out the wall) and 
also on the chemistry specifics of the wall and impuri- 
ties. Let us focus on the mutual interplays between n and 
z e and in obtaining an equation for their evolution with 
time as flow passes through the channel. In particular, the 
interdependence z e (n) may be naturally expected to be 
continuously decreasing when n increases, to take a finite 
value zo at n = 0 (clean filter), and to saturate to zero 
when n reaches some critical value n sat at which all active 
centers of the wall become well covered by impurities. 
We thus postulate the simplest z e (n) dependence fulfilling 
such conditions: 

7 = if«<« sat 

where the notation || . . . || stands for min{l, . . .}. Obvi- 
ously, other sensible choices for z e (n) are possible such as, 
e.g., z e (ri) = Zo(l — \\n/n sat \\) Zl , withzi as a positive coef- 
ficient that probably covers at a good approximation most 
actual possibilities depending on its value, as it corre- 
sponds to a z e (n) functionality with downwards curvature 
if 0 < z\ < 1, to no curvature for z\ = 1 (i.e., Equation 1), 
and to upwards curvature for z\ > 1. For simplicity, we 
shall consider in this letter only the case z\ = l. Note also 
that z e may depend on position and time via the n(x, t) 
dependence. 

Impurity trapping probabilities as a function of z e and n 

The role played by z e in our model will be in fact twofold. 
First, it affects how large the distance is within which if 
the impurity approximates the inner wall then the latter 
attracts the former so much as to consider it as a collision. 
This attraction distance may be seen as an effective radius, 
p e , of the impurity (see Figure 1), so if the distance from 
the center of the impurity to the center of the channel is 
larger than r e — p e , the impurity will actually touch the wall 
(dressed with already trapped impurities). Let us discuss 
the p e (z e ) dependence. We consider first the simplest case 
of an unscreened electrostatic interaction, in which the 
potential energy of an impurity at a distance p e from the 
wall is oc z e /p e . Its kinetic energy associated to the ther- 
mal agitation is oc I<b T. By equating both and also taking 
into account the finite bare size of impurities, we obtain 



p e = p 0 + p^z e as a reasonable approximation, where p[ 
is a constant inversely proportional to temperature. More 
interesting is the case in which ions in the carrying fluid 
partly screen out the electrostatic interaction. The precise 
algebraic distance dependence of the screened electro- 
static energy may be different for each specific channel, 
fluid, and impurity, but we adopt here the common Debye- 
Hiickel approximation in which this energy at a distance 
p e from the surface is taken as oc (z e /p e ) exp(— p e Ad) 
where Ad is the so-called Debye length. In aqueous liquids, 
Ad is a function of the ionic strength, and for concreteness, 
we will consider it to be dominated by the background 
electrolytes in the fluid rather than by the impurities to be 
filtered out (this seems to be the case at least of the mea- 
surements in [5,6]), so is essentially independent on the 
concentration of the impurities to be trapped by the chan- 
nel walls. By equating now the screened potential energy 
at p e to the thermal kinetic energy, we get 




In the right-hand side of this equation, for convenience, 
we have expressed the thermal kinetic energy in units of 
the unscreened potential in the clean channel at a distance 
po from the surface, so p\ is an nondimensional coefficient 
proportional to T. We have also taken into account the 
finite bare size of the impurities by using p e — po instead of 
Pe in the potential energy term. From the above equation, 
Pe can be obtained with the help of the principal Lambert 
W function as follows: 

Pe = P0 + XDW (^\. (3 ) 

Although W(x) can be easily evaluated by modern com- 
puters, it is worthwhile to mention its asymptotes W(x) ~ 
x for x 1 and W(x) ~ Inx — lnln# for x > 3. 
In particular, the first limit means that for small values 
of z e , it is recovered in the linearity of p e with z e found 
for the unscreened interaction, with p[ = po/(piZo). In 
the remaining part of this letter, we shall use the full 
Equation 3 for the p e (z e ) functionality. 

We may now obtain the fraction f e of impurities that 
flow, at given t and x values, near a collision distance 
from the impurity-dressed wall. For that, we assume 
that the fluid velocity profile is given by the Poiseuille 
law, [10] u{r) oc rl — r 2 , where u is the fluid veloc- 
ity and r the distance to the channels axis (see [11] for 
an explicit discussion supporting that at least for chan- 
nels of radius > 10 nm, the flows of water-like liquids 
driven by hydrostatic pressure are in fact in the Poiseuille 
regime). Then,/ e is given by the fraction of the fluid mass 
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that passes through the outer ring r e — p e 



< r < r e 



i.e.,f e = fr*- Pe u ( r ) r dr/ u(r)rdr. The result of those 
integrations is 



(4) 



In the considerations leading to Equation 4, we have 
implicitly taken the concentration of impurities as con- 
stant along the radial coordinate r. However, in principle, 
it could be expected that near the walls the electric poten- 
tial will influence the distance between impurities. To test 
whether this effect may be of relevance, a Debye-like con- 
centration profile C imp (r) = C bu i k (l + z e C / £T (r " re)Az >) 
was also considered. The corresponding f e is then given 
by f^_ pe Q mp (r)u(r)rdr/ Q mp (r)u(r)rdr, the explicit 
algebraic result being too cumbersome to be reproduced 
here. As it will be commented on in detail later in this let- 
ter, we have observed that both Equation 4 and the more 
complicated alternative are able to predict essentially the 
same filtering performances and time evolutions, and so 
in the following, we will employ the simpler Equation 4 
unless stated otherwise. 

The second influence played by z e in our model con- 
cerns the probability that an impurity gets actually bound 
to the inner wall of the channel once it actually is within 
a collision distance from that wall. We express the proba- 
bility that a given impurity entering a differential slice of 
the channel with thickness dx gets trapped in that slice as 
^trap = /entrap dx> where £^ rap is then a trapping prob- 
ability per unit length for the impurities flowing near a 
collision distance from the surface. This £^ rap will obvi- 
ously depend on the chemistry of impurities and active 
centers of the nanostructure and also on the number den- 
sity of active centers not yet saturated by existing bindings. 
The latter indicates that £^ rap will grow with z e , and in 
particular, we may adopt the natural first-order approx- 
imation £2[ rap = + &iz e (^o corresponds then to 
the value in a conventional non-nanostructured filter and 

Equation for Bn(x, t)/Bt 

Let us now build, on the basis of the above relationships, 
equations for the evolution of the areal density of trapped 
impurities, n, as a function of time t and position x when 
an impure fluid flows through the channel due to hydro- 
static pressure. We start by considering the differential 
channel slice going from x up to x + dx, and we write down 
the expression for dn/dt evaluated at that slice, which 
may evidently be expressed as Oi mp £2 t rap/(27rrnd#), where 
^imp is the flow of impurities brought by the incoming 
fluid (in units of s _1 ; the factor (27rrnd#) _1 is due to the 
areal density normalization in the definition of n). In its 
turn, Oi mp can be written as Oi mp = Q mp O where O is 



the fluid flow and Q mp the incoming number concentra- 
tion of impurities. Gathering together the previous results 
in this letter, we get 



C imp (*,QO(0 r , 0 -I 



p e (z e (n(x, t))) 



r 0 - nn(x, t) 



(5) 



with the z e (n) and p e (z e ) dependences given by 
Equations 1 and 3. 

Equations for <!>(£) and &Q mp (x, t)/Bx 
In order to solve the filtration dynamics (i.e., to obtain 
n(x,t) and Q mp (#,£))> it is necessary to supplement 
Equation 5 with formulas for 0(£) and Q mp (#, t). Regard- 
ing the fluid flow, we apply the Poiseuille law for incom- 
pressible fluids of viscosity r] in a cylindrical channel of 
length L and radius r e (x, t): [10] 



0(£) 



nP 



8r) j Q [ro — rxnix, t)\ dx 



(6) 



In this equation, P is the pressure difference between 
both ends of the finite-length channel, which we take con- 
stant with time. Note that O becomes zero when at some 
x, the n value becomes n clog = rn/ri, i.e., r e becomes 
zero at that location and the channel becomes fully closed 
by impurities. Note also that Equation 6 reduces in the 
particular case ro ^> r\n(x, t) (which is common in 
experiments) to O = 7TPrQ/(8rjL). 

We construct now the supplementary equation for 
C[ m p(x,t). For that, we again consider the differential 
channel slice going from x up to x + dx. The number 
of impurities that become trapped in its walls during 
an interval dt is (27zrodx)(dn/dt)dt (the factor Inrodx is 
again due to the areal normalization in the definition of n). 
The numbers of impurities entering and exiting the slice 
in the liquid flow are 0(£)Q mp (#, t)dt and 0(£)Q mp (# + 
dx, t)dt respectively. Mass conservation balance therefore 
gives 



^ ^imp (Xf t) 
dx 



-27rro dn(x, t) 



oa) 



dt 



(7) 



Notice that Equations 5 to 7 are coupled to each other. 
In fact, they form now a closed set that can be numerically 
integrated by providing the specific values for the charac- 
teristics of the filter, for any given pressure difference P 
and incoming impurity concentration Q mp (0, t). In what 
follows, for simplicity, we will always consider for the lat- 
ter a constant value Co. The computation to numerically 
integrate Equations 5 to 7 is relatively lightweight (e.g., 
calculating our Figure 2 took about 15 min in a current 
personal computer that considered 2 x 10 4 finite-element 
^-slices). 
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Figure 2 Time dependence, (a) Results, obtained by integrating Equations 5 to 7, for the time dependence of the areal density of trapped impurities 
(continuous lines) at the entrance of the channel n(x = 0, t) and at its exit point n(x = L, f), and also the global average areal density of trapped 
impurities n(f) = L~ ] n(x, t)dx. The areal density axis is normalized by the saturation value n sat . The time axis is normalized by the half-saturation 
time, defined by n(ti / 2 ) = n sat /2. The parameter values used are as follows (see main text for details): p 0 = 1 3 nm, pi = 0.1 1, ko = 5.1 nm, 
n sat = 1.5 x 10 15 /m 2 = 0.25/ p^r 0 = 500 nm,ri = 8.8 x 1(T 24 ,m 3 = 4p^Q 0 = 0, Q]Z 0 = 1.2 x 10 5 /m,/. = 7.25 mm,P = 3 x 10 5 Pa, r\ = 10 -3 
Pa-s, and Cj mp (x = 0, f) = 1 0 10 /m 3 . We also show the linear, logarithmic, and saturated behaviors (as dashed, dotted, and dot-dashed lines 
respectively), (b) Time dependence of the logarithmic removal value (LRV), calculated using the same parameter values as in Figure 2a. 



Discussion of the results obtained by integrating 
the model equations 

Numerical integration and comparison with some existing 
partial measurements 

We show in Figure 2 an example of the results obtained by 
numerically integrating Equations 5 to 7 using some rep- 
resentative values for the parameters involved (and always 
in the case of constant P and Q mp , and starting from a 
clean initial state n(x,t = 0) = 0). In particular, we have 
chosen parameter values that reproduce the case of chan- 
nels coated with Y2O3 nanopowders as measured in [5] 
(they are essentially valid also for the quite similar case 
of channels with ZrC>2 nanocoating reported by the same 
group in [6]). In these niters, the channels have a typical 
value of the nominal radius ro = 500 nm and length L = 
7.25 mm. They were shown [5] to efficaciously retain MS2 
viruses (of radius po = 13 nm) carried by water with NaCl 
as background electrolyte and a conductivity of 400/xS/cm 
(corresponding then to Ad — 5.1 nm) feed at a pressure 
P = 3 bar. The incoming impurity number concentra- 
tion was Q mp (x = 0,0 = 10 10 /m 3 . For the saturation 
areal density n sat , we will estimate, based on figure nine of 
[5], a quite conservative value n sat = 1.5 x 10 15 /m 2 , cor- 
responding to 0.25/Pq. For the parameter n, we will use 
the value 4pg, also consequent in the order of magnitude 
with figure nine of [5]. These numbers imply that at sat- 
uration (n = n sat ), the effective radius of the channel is 
r sat _ ^gy nm ]\j 0 t e triat j-hjg va i ue i s rather close to the 

clean-state value of 500 nm, and then it would correspond 
to an increase of the hydrodynamic resistance of only 
about 10% (unfortunately, the nanocoatings in [5,6] seem 
to be washed out before they can be fully saturated; how- 
ever, other nanocoated filters [4,7,8] have been shown to 
have hydrodynamic resistance only moderately increased 



at saturation, what is indeed an advantage of paramount 
importance for applications). We will also assume a null 
^trap va l ue at the saturated state, i.e., £2o = 0 (so that we 
neglect conventional filtration mechanisms and focus on 
the effects of nanostructuring alone). In order to proceed 
with the numerical calculation of Equations 5 to 7, only 
two parameters remain to be given estimated values: Q\Zo 
(Qi and zo do not appear separately in Equations 5 to 7) 
and pi (or equivalently, via Equation 3, the effective impu- 
rity radius in the clean state of the channel, Pg lean ). We 
have found that the values ^1^0 = 1.2 x 10 5 /m and p\ = 
0.11 produce results in reasonable agreement with the 
available experimental information, as we discuss below. 
The value p\ = 0.11 corresponds to pf e * n = 33 nm, or 
po + 4A D . 

Figure 2a presents the results corresponding to integrat- 
ing Equations 5 to 7 using these parameter values, for 
the areal density of trapped impurities at the entrance of 
the channel n(x = 0, t) and at its exit point nix — L, t) 
and also for the global average areal density of trapped 
impurities h{t) = L~ l n(x, t)dx. Figure 2b presents the 
corresponding logarithmic removal value (LRV), calcu- 
lated as — log 10 [ Q mp (x = L, t)/Q m p(x = 0, £)]. Note that 
in Figure 2a,b, the time axis is logarithmic and that for 
convenience, it was normalized by the time £1/2 denned by 
the condition #(^1/2) = n sat /2 (half- saturation time). The 
agreement of these numerical results with the measured 
filtration performance reported in [5,6] is fairly good. In 
particular, we obtain an initial LRV of 6.5 log, equal to the 
LRV measured in [5,6] when the actual filters (composed 
by a macroscopic array of microchannels) were chal- 
lenged with only about 1 L of water (the authors of [5,6] 
estimate that such volume carries a total amount of impu- 
rities that is orders of magnitude smaller than the total 
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available binding centers in their filter, so the measure- 
ment is expected to correspond to almost clean channels, 
as in fact seems to be confirmed by microscopy images 
[5]). The calculated LRV is of 4 log at t/ty 2 - 0.7, which 
is also in fair agreement with the observation of a 4 log 
nitration in [5,6] after passing through the macroscopic 
filter approximately from 200 to 1,000 L, depending on 
the measurement. However, obviously, a more stringent 
determination of the parameter values, and in general of 
the degree of validity of our equations, would need more 
precise and detailed data. Unfortunately, to our knowl- 
edge, no measurements exist for the time evolution of the 
filtering efficiency of channels with nanostructured walls 
with a ^-density and precision sufficient for a fully unam- 
biguous quantitative comparison with the corresponding 
results of our equations; in fact, one of the main moti- 
vations of the present Nano Idea Letter is to propose 
(see our conclusions) that such measurements should be 
made, in order to further clarify the mechanism behind 
the enhanced impurity trapping capability of the channels 
with nanostructured inner walls. 

As a further test, we have repeated the same numer- 
ical integration as in Figure 2a,b but considering a 
radial impurity concentration profile Q mp (r) oc 1 + 
z e C'e~( r ~ re ^ kD , instead of a constant one as in Equation 4. 
We have obtained very similar results, provided that the 
parameter Q\Zo is conveniently varied: In particular, we 
observed that the filtration dynamics results obtained 
using Equation 4 and any given value y for Q\Zo can be 
reproduced using the above Debye-like profile if employ- 
ing for Q\Zo a new value (specifically, the new value can 
be estimated, by comparing the initial nitration perfor- 



P (r) 



mance, as y2^_^ dean 2^ mp ^V^o J ^dL» where ^ = 
f^ e y(x)u(x)dx; for instance, taking zo C = 200, which 
probably is a fair first approximation for the measure- 
ments in [5-8], the parameter values used in Figure 2 
correspond to 3.2 x 10 4 /m as equivalent Q\Zo value when 
using the Debye approach). These results indicate then 
that, as it could be expected, the wall charge effect on the 
radial gradient of the concentration may be safely sum- 
marized, for our present purposes, as one of the factors 
influencing the value of the trapping probability coeffi- 
cient Q\Zo to be used when applying Equation 4. 

Linear, logarithmic, and saturated approximations 

In Figure 2a, it is possible to identify in our results for the 
areal density of trapped impurities some ^-ranges in which 
the ^-dependence is relatively simple: (1) The initial time 
behavior is an approximately linear n(t) growth; (2) in the 
intermediate regime, the growth of n(t) becomes approx- 
imately logarithmic; and (3) at sufficiently large t values, 
the saturation limit is reached, in which n approaches a 
value n sat at a slow pace. These regimes are easily seen in 



Figure 2a for n(x = 0, t), n(x = L,t), and h(t), albeit in 
each case they are located at different t/t\/2 ranges. The 
figure also evidences that it is possible for the linear and 
logarithmic ^-ranges to overlap each other (the case of nit) 
with the parameter values used in Figure 2). 

In the case of a very short cylindrical channel (so that 
all ^-derivatives may be neglected), it is possible to find 
analytical expressions for the nit) evolution in the lin- 
ear and logarithmic regions: For the linear regime, by just 
introducing in Equation 5 the condition t ~ 0, we find: 



nit) ~ A Xm t 



with 



^lin = 



(Qp + ^iZo)C imp P 

16r]Lro 

Po~ro + k D W | 



Po 



(8) 



(9) 



The logarithmic regime can be found by using the con- 
dition n ~ n sat /2: 

n sat t 
nit)- — +Ai og ln— - (10) 

2 h/2 

with 



^log — 



(2^o + £2 1^0 ) h/2Ci m pP 



(2p 0 - 2r 0 + l D W (j^j + nn s ^ - (2r 0 - nn Si 



(ID 



In obtaining the above Equations 8 to 11, we have 
assumed that n(0) =0 and that p e < r e at t = 0 or t\/2. 

Conclusions and proposals for future work 

This letter has proposed a model for the main generic 
features of the channels with nanostructured inner walls 
with respect to trapping and accumulation of impurities 
carried by fluids. This includes, e.g., their capability to 
clean the fluid from impurities of a size much smaller 
than the channels' nominal radius, with comparatively 
small resistance to flow (much smaller than in conven- 
tional channels with a radius as small as the impurities). 
The model attributes the enhanced nitration capability to 
the long-range attraction exerted by the exposed charges 
in the nanostructured walls and also to their binding capa- 
bility once the impurities actually collide with them. Both 
features were quantitatively accounted for by means of a 
phenomenological effective -charge density' of the nanos- 
tructured wall. The model also predicts the time evolution 
of the trapped impurity concentration and of the filter- 
ing capability, including three successive regimes: a linear 
regime, a logarithmic regime, and the saturated limit. 

We believe that our equations could make possible some 
valuable future work, of which two specific matters seem 



Ramallo Nanoscale Research Letters 201 3, 8:1 9 
http://www.nanoscalereslett.eom/content/8/1/19 



Page 7 of 7 



to us more compelling: First, it would be interesting to 
check at the quantitative level the agreement with exper- 
iments of the time evolutions predicted above. For that, 
we propose to perform time-dependent measurements 
made in controlled flow setups. We have chosen in our 
equations flow constraints which seem appropriate for 
this purpose; however, the model can be tested with dif- 
ferent setups by just numerically integrating Equation 5 
under the corresponding experimental constraints. One of 
the prime purposes of this letter is in fact to suggest such 
measurements. 

A second interesting future work (already in progress in 
our research group) is the design of optimal geometries for 
the combinations of channels forming filters. Our model 
opens this possibility because of the explicit use of the 
Poiseuille relations that allow the calculation of the resis- 
tance to flow of complex associations of those channels, 
in series and/or parallel. The effective diffusivity and tor- 
tuosity of the pathways' network are also accounted for by 
these equivalent-circuit analyses. 
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